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Abstract 

Aims. A growing number of previously hidden Galactic X-ray sources are now detected with recent surveys performed by the 
INTEGRAL and Swift satellites. Most of these new sources eluded past surveys due to their large local X-ray extinction and con- 
sequent low soft X-ray flux. The Swift-B AT performs daily monitoring of the sky in an energy band (15-150 keV) which is only 
marginally affected by X-ray extinction, thus allowing for the search of long periodicities in the light curve and identification of the 
nature of the X-ray sources. 

Metliods. We performed a period search using the folding technique in the Swift-BPS light curves of three INTEGRAL sources: 
IGRJ05007-7047, IGRJ13 186-6257 and IGR J17354-3255. Their periodograms show significant peaks at 30.77 ± 0.01 d, 
19.99 ± 0.01 d and 8.448 ± 0.002 d, respectively. We estimate the significance of these features from the x' distribution of all 
the trials, finding a probability < 1.5x10""* that the detections occurred due to chance. We complement our analysis with the study of 
their broadband X-ray emission. 

Results. We identify the periodicities with the orbital periods of the sources. The periods are typical for the wind accretors X-ray 
binaries and we support this identification showing that also their energy spectra are compatible with an X-ray spectral emission 
characteristic of high-mass X-ray binaries. The spectrum of IGR J05007-704 that resides in the Large Magellanic Cloud, does not 
show any intrinsic local absorption, whereas the spectra of the Galactic sources IGR J17354-3255 and IGR J13186-6257 may be 
affected by a local absorber The folded light curve for IGR J 13 186-6257 suggests a possible Be companion star. 
Conclusions. 
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1. Introduction 



X-rays: individual: IGR J05007-7047 - X-rays: individual: IGR J13186-6257 - X-rays: individual: 



One 



of the main results of the IBIS/ISGRI telescope 
20031) on board the INTEGRAL satellite 



( Ubertini et al 



(iWinkler et al . 2003) is the detection of a large number of new 
hard X-ray sources. At present, INTEGRAL has detected more 
than 500 sources (~ 26% are associated with extragalactic 
sources, ~ 26% with Galactic sources and ~ 48% whose nature 
has not yet been determined!]. Among the Galactic sources, 56 
are identified as binary systems with a supergiant companion. 
Most of them, characterized by persistent emission, are likely 
deeply embedded in their stellar wind, as sugges ted by their 
strong intrinsic absorption (Nh > 10^^ cm^, see e.g. IWalter et al.l 
12001 



An important contribution to the study of these new 
Galactic sources is pro vided by the Burst Alert Telescope (BAT , 
iBartheknv et"an l2005h on board Swift dGehrels et al.l l2004 . 
which has been performing a continuous monitoring of the sky 
in the hard X-ray energy range (15-150 keV) since November 
2004. The telescope, thanks to its large field of view (1.4 stera- 
dians half coded) and its pointing strategy, covers a fraction be- 
tween 50% and 80% of the sky every day. This has allowed 
the detection of many of the new INTEGRAL HMXBs (e.g. 
[Cusumano et al. 2010b) and the collection of their long term 



' An updated list of these sources with their main 
properties and the relevant references can be found in 
[http ://irfu.cea.fr/Sap/ IGR-Sources/| 



light curves and of their averaged X-ray spectra. The long and 
continuous monitoring of these sources allows to investigate the 
intrinsic emission variability, to search for long periodicities (or- 
bital periods) and to discover the presence of eclipse events. The 
role of ^wi/f-BAT is therefore fundamental to unveil the nature 
and the geometry of these X-ray binary systems. 

We are performing a systematic study of the light curves of 
the new INTEGRAL sgHMXB sources to search for > 0.5 day 
periodicities and to discover the presence of eclipse events al- 
lowing the study of the binary system geometry. In this paper 
we present the detailed timing analysis and broadband (0.2-150 
keV) spectral analysis of IGR J05007-7047, IGR J13 186-6257, 
and IGRJ17354-3255. 

IGR J05007 -7047 (also known as IGR J05009-7044), lo- 
cated in the Large Magellanic Clou d (LMC), was d etected with 
INTEGRAL in the 17-60 keV band dSazonov et al.ll2005.) with a 
flux of 1.2 X 10~" erg cm"^ s ' , corresponding to a luminosity 
of 3.6x10"^^ erg s"', assuming a distance of 50kpc ("Guin an et alJ 
1998). A follow-up Chandra observation allowed its association 
with the relatively bright blue (y =14.8 mag, B - y=-0.01 mag) 
star US NO-B 1.0 01 92-005 7570 (lHalpernl2005l) of spectral type 
B2 III (iMasetti et alj|2006l) at a redshift consistent with that of 
LMC. iLa Parola et al.l ( 2010 b) reported on the detection in the 
BAT surv ey data of a perio dic modulation at 30.77 d, later con- 
firmed by ICoe et all (1201 Ol) through the analysis of the optical 
counterpart light curve. 
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IGR J 1 3 1 86-6257 was fi rst reported in the Third Integral 
Catalogue (iBird et alJ 12007 *) as an unidentified hard X-ray 
source. A follow-up Swift-XRT observation better constrained 
its p osition and revealed a flat soft X-ray spectrum (iLandi et al.l 
12008.) . A short (~ 5 ks) Chandra observation found an X-ray 
source wifliin the Swift-XRT error box at (RA, Dec [12 0001= 
199.604500 deg, -62.970972 deg) (Tom sick etan i2009). The 
Chandra source, CXOU J131825. 0-625815, was identified as 
the most likely counterpart of the INTEGRAL source, with a 
low spurious detection probability of ~ 5%. The Chandra po- 
sition allowed for the counterpart optical identification with the 
2MASS source J13182505-6258156. The spectral shape of the 
source was consistent with a highly absorbed (1.8 x lO^-' cm"^ ) 
power-law spectrum. 

IGRJ17354-3255 was discovered during the INTEGRAL 
Galac tic Bulge monitoring program (Kuul kers et al.l l2006i 
120071) . The source has also a possible, highly vari- 
able, y-ray counterpart detected with the AGILE satellite 



(AGLJ1734-3310, iBukarelli et al.1 
follow-up observation 



20091). A Swift-XRT 
(IVercellone et al .1 120091) revealed the 
presence of two sources within the INTEGRAL error circle. 
Based on its variability, the authors suggested the source at (RA, 
Dec [J2000]= 263.863167 deg, -32.930250 deg) as the likely 
counterpart of IGR J 17354-3255. The spectrum of this source 



is characterized by strong absorption (Nh =5+2 x 10 cm" ) 
and high flux variability (an order of magnitude change in flux 
within 7 hours). A further Chandra observation dTomsick et al.l 
12009 ) confirmed its variable nature and suggested the associa- 
tion with the IR source 2MASS J17352760-3255544. A flux 
modulation in the Swift- BAT data, w ith a period of 8.452 days, 
has been first reported bv lD'Ai et al1 C2010). 

This paper is organized as follows. Section 2 describes the 
data reduction. In Sect. 3 and 4 we describe the temporal and 
spectral analysis and in Sect. 5 we discuss our results. 



2. Data Reduction 

The raw BAT survey data of the first 54 months of the Swift 
mission were retrieved from the HEAS ARC pubhc archivaj 
and processed with a dedicated software (ISeereto et al.l 1201 OIK 
that performs screening, mosaicking and source detection on 
^wi/f-BAT data and produces spectra and light curves for any 
given sky position. IGR J05007-7047, IGR J13 186-6257 and 
IGRJ17354-3255 were detected in the Swift-BAT all-sky map 
with a significance maximized in the 15-50 keV band of 16.1, 
10.0, and 18.9 standard deviations, respectively. Light curves 
were extracted in the 15-150 keV energy range with the maxi- 
mum available time resolution (~ 300 s) and corrected to the so- 
lar system barycentre (SSB) by using the task EARTH2SUN. The 
15-150 keV spectra were obtained extracting the source fluxes 
from the sky maps in sixteen energy bands, and were analyzed 
using an appropriately rebinned CALDB response matrix. 

Swift-XRT observed IGR J05007-7047 between 2010-10- 
25 and 2010-10-28 (Obsid 31846); IGR J13 186-6257 was 
observed on 2008-05-01 and on 2010-05-25 (ObsId 37071); 
IGRJ17354-3255 was observed on 2008-03-11 and 2009-04- 
17 (ObsId 37054). The sources were all observed in Photon 
Counting mode. Source and background spectra and the relevant 
ancillary files were created and retrieved from the Swift-XRT 
product on-line builder mai ntained by the Sw ift data center at 
the University of Leicester (lEvanset al .1120091) . The Swift-XRT 
spectra were rebinned to have at least 20 counts per bin, and 



Swift-BAJ" energy channels were grouped to obtain a S/N ratio > 
3 per channel. This choice allows us the use of the;^'^ statistics. 
In the following, we report errors at 90 % confidence level, if not 
stated otherwise. 

Table [Urepoits the log of the Swift-BAT and Swift-XRT ob- 
servations used for the present analysis. 

3. Temporal Analysis 

We analyzed the long term ^vvZ/f-BAT light curve to search 
for intensity modulations by applying a folding technique and 
searching in the 0.5-lOOd period range. The period resolu- 
tion is given by P^/(A^ ATbat), where = 16 is the number 
of trial profile phase bin s and ATbat is the data span length 
(iBuccheri & Saccolll985l) . 

The average rate in each phase bin was evaluated by weight- 
ing the light curve rates by the inverse square of the correspond- 
ing statistical error 



Ri 



I. I /erf 



(1) 



|http : / /heasarc .gsfc.nasa. gov/ docs/ archive . html 



where Rj is the average rate in the j-th phase bin of the trial 
profile, rj are the rate of the light curve whose phase fall into 
the j-th phase bin and er,- are the corresponding statistical errors. 

The error on Rj is ^ ^JYAJerj^ . The weighting procedure was 

adopted to deal with the large span of er, and it is justified be- 
cause the BAT data are background dominated. 

The significance of a feature P„ in the resulting periodogram 
cannot be evaluated using the statistics because, if the source 
is variable on time scales comparable to those under investiga- 
tion, the average is far from the value expected for white noise 
N - 1). Therefore, the significance of Pq is eva luated as follows 
iLa Parola etaDbOlOaLlCusumano et al.ll2010al) : 

1 . We fit each periodogram with a second order polynomial and 
subtract the trend to the distribution, obtaining a new xf: dis- 
tribution. 

2. We build the histogram of the new;^'^ distribution excluding 
the interval around P,, and those around any features deriving 
from it (e.g. periods multiple of Po). 

3. We fit the tail of the resulting distribution with an exponential 
function and evaluate the integral of the best-fit function beyond 
Xci^o)- This integral yields the number of chance occurrences 
due to noise. 

3.1. IGRJ05007-7047 

Figure [T^ shows the periodogram obtained for 
IGR J05007-7047. We find significant evidence for the presence 
of a periodicity 0^^^ ~ 210) at a period of Pq = 30.77 ± 0.03 d, 
where the period uncertainty is the period resolution at Pq. The 
other two significant features in the periodogram correspond 
to 2Po and 3Po. The value of Xc (i-^- ^ft^'" subtracting 
the periodogram trend) at P^ is 177. Figure [TJi shows the Xc 
distribution excluding the values around P^ and its multiples. 
We then fitted the tail (xf, > 15) of the resulting distribution with 
an exponential function and evaluated the integral of the best-fit 
function beyond P,,. The number of chance occurrences to have 
a value of 177 or higher is 3.9 x 10"'^, that corresponds to ~ 7 
standard deviations in Gaussian statistics. The pulsed profile 
(Fig.[l];) folded at P„ with Tepoch = 54159.753 MJD, shows a 
roughly sinusoidal modulation with a minimum value consistent 
with zero intensity. The centroid of the minimum, evaluated 
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Table 1. Log of 5wj^-BAT and 5w//r-XRT observations used for this work. 



by fitting the data around the dip with a Gaussian model is 
at phase 0.23. The uncertainty depends on the rough binning 
of the folded light curve and we conservatively assume an 
error corresponding to one bin phase; the dip phase, therefore, 
corresponds to MJD (54166.8 + 1.0) + nxPo MJD. The phase 
coverage of the light curve bins shows that this dip is not due 
to an accidental under-sampling of the light curve at this phase. 
We self-consistently checked this aspect also for the dips of 
IGR J17354-3255 and IGRJ13 186-6257. 

The three Swi/f-XRT observations were performed at the or- 
bital phase intervals of 0.369-0.380, 0.428-0.432 and 0.469- 
0.482, respectively. The 0.2-10 keV light curve shows a per- 
sistent emission with a variability of a factor ~ 5 and an aver- 
age count rate of 0.058+0.004 count s"'. The source events se- 
lected in the three observations are 77, 37 and 509, respectively. 
Because of the low statistic content of the first 2 observations we 
restricted the analysis to Obsld 00031846003. 



3.2. IGR J1 31 86-6257 

The highest feature, with a value of ~ 164, is at - 
19.99 + 0.01 d (Fig. |2ji) . The second and third highest peaks 
in this periodogram are found at multiple values of P„ at ~ 40 
d and 60 d. We fit the periodogram with a second order polyno- 
mial and subtracted the trend from the distribution, finding a 
X^=126 at ?„. We therefore built the histogram of the Xc distri- 
bution (Fig.|2]3) excluding the interval around P„ and the first two 
multiples of Pq. The integral above Pq yields a number of chance 
occurrences due to noise of 3.2 x 10 corresponding to a sig- 
nificance for the detected feature of ~ 4.7 standard deviations in 
Gaussian statistics. 

The pulsed profile (Fig. |2};) folded at Pg with 
Tepoch=54417.993 MJD, shows a plateau between phase 
0.9 and 1.4, where source emission is dominated by back- 
ground, whereas the source becomes clearly detected in the 
remaining half of the orbit. We evaluated the time passage 
corresponding to the folded light curve minimum bin at phase 
0.95 con-esponding to Tdip = (54436.2 + 0.6) + nxP^ MJD. The 
two Swift-XRT observations were performed at orbital phases 
0.50 and at phase 0.20. The source is clearly detected in the first 
observation with an averaged count rate of 0.28+0.02 count s"' . 
In the second observation, performed two years later, no source 
is detected at the same celestial coordinates. We estimated an 
upper limit on the source flux between 1.2 and 1.5 x 10"'^ erg 
cm"^ s ' depending on the choice of the continuum as derived 
from the first observation. For spectral analysis we made use 
only of the data from the first observation. 



3.3. IGR J1 7354-3255 

The periodogram obtained for IGR J17354-3255 (Figure [3^) 
shows significant evidence for the presence of a periodicity 
(X^ ~ 121) at a period of Pq = 8.448 ± 0.002 d. The other 
significant features that appear in the periodogram corresponds 
to multiple of Po (up to 8Po). After correcting for the pe- 
riodogram trend, we obtain ;t'^(Po)=94. Figure [3j3 shows the 
X^ distribution excluding the values around P,, and its multi- 
ples. The number of chance occurrences to have a value of 
94 or higher is 1.5 x 10""* that corresponds to ~ 4 standard 
deviations in Gaussian statistics. The pulsed profile (Fig. [3};) 
folded at P,, with Tepoch=54175.159 MJD, shows a modulation 
with a minimum value consistent with zero intensity. The cen- 
troid of the minimum is at phase 0.78 corresponding to MJD 
(54181.75 ± 0.26) + nxP,, MJD. 

The two Swift-XRT observations were performed at the or- 
bital phase intervals of 0.719-0.767 and 0.304-0.338, respec- 
tively. The phase of first observation corresponds to the mini- 
mum in the folded profile (Figure|3];), while the phase of the sec- 
ond one falls around the maximum of the profile. Figure|4]shows 
the Swift-XRT field of view of the two observations: th e posi- 
tion of the two possible counterparts reported in Vercellon e et al.l 
(l2009h are marked as Srcl (RA, Dec [J2000]= 263.863167 deg, 
-32.930250 deg) and Src2 (RA, Dec [J2000]= 263.826333 deg, 
-32.906361 deg). Srcl is not detected in the first observation, 
with a 3 cr upper limit of 2.3 x lO""* count s"'. As this observa- 
tion corresponds to a phase consistent with the minimum of the 
BAT folded light curve we can confidently associate Srcl to IGR 
J 17354-3255. In the second observation the 0.2-10 keV average 
count rate of Srcl is 0.076+ 0.004 count s"', with a variability 
of a factor ~ 20. 

4. Spectral Analysis 

We have performed a broad band spectral analysis of the XRT 
and BAT data. The XRT and BAT spectra are not simultaneous, 
as the BAT spectra are accumulated over 54 months. Therefore, 
we have checked for the presence of spectral variability in the 
BAT data by building the hardness ratio defined as the ratio be- 
tween the count rate in the 15-150 keV and in the 15-35 keV 
energy band, with a time resolution of 1 day. We found no sig- 
nificant variation of the hardness ratio along time for the three 
sources. We also verified that the BAT count rate during the XRT 
observations is consistent with the average count rate in the 54 
months. However, to take into account residual flux and inter- 
calibration uncertainties between the two spectra, we included a 
multiplicative constant factor (Cbat) in all the models. We kept 
it frozen to unity for XRT data and free to vary for BAT data. 
Because of possible correlations between the values of Nh and 
the indeces of the power-laws, we built 2D contour-plots around 
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the best-fitting values, always checking the consistency of the er- 
ror estimates. We report in Tabj2]the spectral best-fitting results 
for the three sources, using a cut-off power-law model. Swift- 
BAT fluxes include the Cbat factor in the calculation. 

4.1. IGRJ05007-7047 

The combined XRT-BAT spectrum was first modelled with a 
simple absorbed power law. The resulting (140.2 with 38 
degrees of freedom, d.o.f.) is not acceptable, and the residu- 
als suggest the presence of a cut-off at high energy. We have 
then included such cut-off (model CUTOFFPL), obtaining an im- 
proved of 32.0 for 37 d.o.f., with a cut-off at Ec = 12 ± 2 
keV and a photon-index r=0.3+0.2. The interstellar absorption 
is consistent with the expected value in the direction of LMC 
(lKalberlaetalJl2Q05h . Because of the unambiguous identifica- 
tion of the optical counterpart as belonging to the LMC, we de- 
rive an extrapolated 0.2-150 keV X-ray luminosity of ~ 9 x 10^'' 
erg s"'. Figure|5^ shows the spectrum, best fit model and resid- 
uals in units of standard deviations. 



4.2. IGRJ1 31 86-6257 

An absorbed power-law model can satisfactorily model the data 
Oif^/d.o.f. =24/21). For this model the Nh is 17+^ x lO^^ cm^^ , 
the power-law photon index is 2.3+0.3. The costant of normal- 
ization between the Swift-XRT and ^wi/f-BAT data is, however, 
inconsistent with unity (being 0.2^Q Qg). We found only marginal 
evidence for a high-energy cut-off in the spectrum (Ec„, > 10 
keV), with a new ;t'^/d.o.f. of 21/20. The position of the high- 
energy cut-off strongly depends on the value of the Nh and of 
the photon-index. Choosing a reference value for Ecu, at 16 keV, 
we found an upper limit to the Nh at 11 x 10^^ cm"^ and a 
photon-index constrained in the 0-1.0 range. Using this model 
the Cbat value is very low (see Tab.|2]i , and this suggests strong 
variability, probably related to the orbital phase. Taking into ac- 
count the time of the Swift-XRT observation, occured at phase 
0.5, we extracted a phase-selected 5w/ff-BAT spectrum, choos- 
ing the interval 0.5-0.8 as representative of the phase where 
source emission is clearly detected above background. In this 
case, the Cbat parameter assumes the value 0.28;^y j2, while the 
spectral parameters are still consistent with those obtained using 
the phase-averaged spectrum. The sensible change in the Cbat 
parameter value supports the idea that the X-ray emission may 
be confined only in a restricted orbital phase interval. We show 
in Fig|5j', the best-fitting model together with data (using the 
Swift-BAT phase-averaged spectrum) and residuals in units of 
cr. 



4.3. IGR J 17354- 3255 

The fit of the combined XRT-BAT spectrum with a simple ab- 
sorbed power-law model gives a statistically acceptable result 
O^^^/d.o.f. =29.6/27), with a Nh value of \Q.6^_\l x 10^^ cm^^ 
and a photon-index of 2.6+0.3. However, fitting the data with 
a cut-off power law, we obtained a significant improvement of 
the fit, with a final ;t'^/d.o.f. of 23.2/26, and an F-test probability 
of ~ 1% of obtaining this improvement by chance. The cut-off 
energy is however poorly constrained and we obtained only a 
lower limit at ~ 18 keV. The photon index F is between 1.0 and 
2.0, while the Nh is 7±2 x 10^^ cm"^ . The expected Nh value 
in the direction of the source is considerably less (1.2-1.6 x 10^^ 



cm"^ ), and this suggests that a fraction of the X-ray absorbing 
medium could be local to the X-ray source. 

The source position is only at 5° from the Galactic Center 
(GC), and if we assume the distance of the GC as the possible 
distance of the source (8 kpc), we obtain an extrapolated X-ray 
luminosity of ~ 5 x 10^^ erg s"'. We show in Fig|5};, the best- 
fitting model together with data and residuals in units of cr. 

5. Discussion and Conclusions 

In the last cens us of high-mass X-ray binaries in the Galaxy 
jLiu et al . 2006). there were 1 14 HMXBs sources reported, with 
only 50 determined orbital periods. With an increasing number 
of new detect ed X-ray sources, our view of HMXBs is rapidly 
changing (see IChatv et al.ll2010l) . also due to the ability of the 
long term operating surveys to search for periodic signals. 

We are currently undergoing a vast program for the 
search of periodic signals exploiting the ^w i/f-BAT surveys 
jLa Parola et al.ll2010atlCusumano et al.ll2010ah . 

In this work, we have shown results from the analysis of 
the data collected by ^wZ/f-BAT during the first 54 months of 
the Swift mission for three still poorly understood INTEGRAL 
X-ray sources: IGR J05007 -7047, IGR J13 186-6257 and IGR 
J17354-3255. 

IGR J05007-7047 was already identified as a possible wind- 
fed X-r ay systern, than ks to the association of the optical coun- 
terpart (lHalpernll2005b . We detect it in the BAT survey at a 
significance level of ~ 16 standard deviations with an average 
flux of ~ 1.3 X 10"^' erg cm"^ s"' in the 15-50 keV en- 
ergy band. The light curve reveals a periodicity of 30.77 ± 0.01 
d that we interpret as the orbital period of the binary system. 
Given the characteristics o f the companion star (spectral type 
B2 III. iMasetti et alj |200^ we can apply Kepler's third law to 
derive the semi-major axis of the binary system. This is given 
by - P^^ X G(M* + M^)IAn'^, where M* and Mx are the 
masses of the supergiant and compact objec t, respectively. We 
adopt Mx = 1.4 M0 and 7 <M* < 12 M© (iHabets & Heintzd 
1981, for a B2 star of radius 4 R© <R* < 5.6 R©). This yields 
84R0 <fl< 98Ro (IIR* <a< 21/?*). The folded profile shows 
a large dip at time of (54166.8 + 0.6) ± nPo MJD with a count 
rate consistent with no emission. However, the width of this dip 
(~ 20 per cent of P,,) is not well constrained, constituting a rough 
upper limit too large to be consistent with the duration of the 
eclipse expected for a circular orbit with the inferred semi -major 
axis. 

The spectrum of IGR J05007-7047 is weU modeUed with 
cut-off power law (F ~ 0.3 and E^. ~ 12 keV). The upper limit on 
the absorbing column is consistent with the Galactic value along 
the line of sight. Therefore this source cannot be classified as an 
absorbed high-mass X-ray binary (> 10^^ cm "^). 

The spectral shape of IGR J13 186-6257 could not be firmly 
constrained, because of the weak statistics of the Swift-XKY 
spectrum. We could not determine any clear cut-off in the en- 
ergy spectrum, with a lower limit at ~ 10 keV. The periodic- 
ity of 19.99 d, is typical for wind/Be accreting X-ray systems. 
The folded light curve shows clear source emission above back- 
ground only for ~ half of the orbital period. Because of the 
poor constraints on the spectral shape, we are not able to uni- 
vocally classify the source. However, the shape of the folded 
light curve, the very low value of the intercalibration constant 
between the pointed Swift-XKI spectrum and the long-term 
^wZ/f-BAT spectrum and the lack of source detection at or- 
bital phase 0.2 suggest a possible Be companion in an eccen- 
tric orbit, where accretion could take place only near the peri- 
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Table 2. Summary of spectral fits using the CUTOFFPL model for the combined Swift-XRT and Swift-BA^ datasets. We report 
absorbed fluxes for the characteristic Swift-XRT (0.2-10 keV) and ^vvZ/f-BAT (10-100 keV) energy bands. E^.,,, parameter for 
IGR J13 186-6257 frozen during the fitting. 



astron passage. The expected equivalent column density value 
in the IGR Jl 3 186-6 257 direction hes betw een 1.24x10^^ cm~^ 
(iKalberla et alJl2005h and 1.59x10^^ cm ^ dPickev & LockmanI 
Il990^ . These estimates are considerably less than our best- fitting 
values, both the simple power-law and for the cut-off power- 
law model. IGR J13 186-6257 may be therefore embedded in a 
denser absorbing medium during accretion phases. Further ob- 
servations are, however, necessary to confirm this scenario . 

IGRJ17354-3255 is detected in the 15-50 keV BAT map 
at a significance level of ~ 19 standard deviations. The average 
flux in this band is ~ 1.56 x 10"" erg cm"^ s"'. The timing anal- 
ysis on the 15-50 BAT light curve unveils an orbital period of 
8.448 + 0.002 d, with a significance of ~ 4 standard deviations, 
confirming its nature as a binary system. The detection of the 
periodicity in enforced by the non-detection of the soft X-ray 
counterpart in the Swift-XRT observation performed at an or- 
bital phase consistent with the minimum of the modulation. The 
folded profile shows a minimum consistent with zero intensity. 
The energy spectrum of IGR J17354-3255 is well modelled by 
an absorbed power law, with Nh ~ 7 x 10^^ cm"^ with a cut- 
off at > 20 keV. This large local absorption suggests that X-ray 
emission from the compact object may interact with the stellar 
wind of the companion star, whose spectral type has not yet been 
determined. 
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Figure 1. (a): Periodogram of Swift-BAT (15-50keV) data for IGR J05007-7047. (b): Distribution of xl values excluding the 
intervals around and its multiples. The continuous Une is the best fit obtained with an exponential model applied to the tail of the 
distribution (x^ > 15). (c): Swift-BAT Light curve folded at a period P = 30.77 ± 0.01 d, with 16 phase bins. 
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Figure!, (a): Periodogram of Swift-BAT (15-50keV) data for IGR J13186-6257. (b): Distribution of xl values excluding the 
intervals around and its multiples. The continuous Une is the best fit obtained with an exponential model applied to the tail of the 
distribution (x^ > 15). (c): Swift-BAT Light curve folded at a period P = 19.99 ± 0.01 d, with 16 phase bins. 



8 



A. D' Ai et al.: The orbital period of three HMXBs 





(c) 





0.5 



1 

Phase 



1.5 



Figures, (a): Periodogram of Swift-BAT (15-50keV) data for IGR J17354-3255. (b): Distribution of xl values excluding the 
intervals around and its multiples. The continuous Une is the best fit obtained with an exponential model applied to the tail of the 
distribution (x^ > 15). (c): Swift-BAT Light curve folded at a period P = 8.452 ± 0.002 d, with 16 phase bins. 
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Figure 4. Field of view of the two XRT observations of IGR J17354-3255, with the position of the two candidate counterparts 
(.Vercellone et al...2009.) . 
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Figure 5. Combined Swift-XRT and Swift-BALT spectra, best-fitting model (from TabjSJi and residuals in units of standard deviations 
of IGR J05007-7047 (upper panel), IGR J13 186-6257 (middle panel) and IGRJ17354-3255 (bottom panel). 



